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Continuum limit of lattice models with Laughlin-like 
ground states containing quasiholes 
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There has been a significant interest in the last years in finding fractional quantum Hall physics 
in lattice models, but it is not always clear how these models connect to the corresponding models 
in continuum systems. Here we introduce a family of models that is able to interpolate between a 
recently proposed set of lattice models with Laughlin-like ground states constructed from conformal 
field theory and models with ground states that are practically the usual bosonic/fermionic Laughlin 
states in the continuum. Both the ground state and the Hamiltonian are known analytically, and we 
find that the Hamiltonian in the continuum limit does not coincide with the usual delta interaction 
Hamiltonian for the Laughlin states. We introduce quasiholes into the models and show analytically 
that their braiding properties are as expected if the quasiholes are screened. We demonstrate 
screening numerically for the 1/3 Laughlin model and find that the quasiholes are slightly smaller 
in the continuum than in the lattice. Finally, we compute the effective magnetic field felt by 
the quasiholes and show that it is close to uniform when approaching the continuum limit. The 
techniques presented here to interpolate between the lattice and the continuum can also be applied 
to other fractional quantum Hall states that are constructed from conformal field theory. 

PACS numbers: 05.30.Pr, 73.43.-f, 03.65.Fd, 11.25.Hf 


I. INTRODUCTION 

Quantum many-body systems can display a wealth of 
peculiar phenomena that are interesting from both a fun¬ 
damental and a practical point of view. An important ex¬ 
ample is topological systems with the possibility to have 
emergent particles with unusual properties like fractional 
charge and nontrivial braiding statistics. The fractional 
quantum Hall (FQH) effect plays a central role in this 
context since topological states were realized experimen¬ 
tally in these systems already in the early eightieJ^l and 
because the states appearing in the FQH effect are sup¬ 
posed to be well-described by relatively simple, analyt¬ 
ical wave functions,I^M^ which is a significant advantage 
for theoretical studies. 

There is currently a lot of interest in investigating pos¬ 
sibilities for obtaining and realizing the FQH effect in 
different systems. This research gives a deeper under¬ 
standing of when and how the FQH effect can occur, and 
it can open up doors to investigate different aspects of 
the effect experimentally since different setups allow dif¬ 
ferent properties to be measured. The hope is also to 
find particularly simple ways to realize the effect, which 
would significantly improve the possibilities for utilizing 
the effect in practical applications. 

In the present paper, we are concerned with FQH mod¬ 
els in lattice systems that are obtained by a strategy, 
in which one tries to keep the wave function as close 
as possible to the corresponding wave function in the 
continuum.l31IHl More specifically, one starts from an an¬ 
alytical FQH wave function in the continuum. A corre¬ 
sponding lattice FQH state is then defined by restricting 
the allowed positions of the particles to a set of lattice 
sites, and it may also be desirable to modify the state 
slightly. Finally, the analytical properties of the state 
are used to find a Hamiltonian for which the state is 


the unique ground state. It has turned out that con¬ 
formal field th eory (CFT) is a very helpful tool in this 
respect lUinillHIll 

The Hamiltonians obtained by using this approach are 
typically valid only for a particular lattice filling factor. It 
is, however, interesting to ask if there is a way to directly 
connect the lattice models to the usual continuum mod¬ 
els, since answering this question would give a more com¬ 
plete understanding of the lattice FQH models, would 
make it easier to compare lattice and continuum FQH 
models, and would provide further guidance on the pos¬ 
sibilities for obtaining FQH models. In [T2l and IT^ it has 
been shown how one can interpolate FQH wave functions 
between the lattice and the continuum limit within the 
context of CFT, but also in these papers a Hamiltonian 
was only found in the lattice limit. Interpolations from 
the continuum towards the lattice usiM an approximate 
Hamiltonian have been done as well,^ but in this case 
the overlap between the ground state of the model and 
the (bosonic) Laughlin state is only high up to a moder¬ 
ate lattice filling factor, i.e. in the quasi-continuum limit. 

Here, we propose an approach that allows us to obtain 
both lattice and continuum models with FQH ground 
states and to interpolate between the two limits within 
the same family of models. The models are built using 
CFT and allow us to compare the CFT Hamiltonian to 
the usual delta function interaction Hamiltonian in the 
continuum. We find that the two do not coincide. The 
construction also allows us to add localized quasiholes to 
the models and study their properties for different lattice 
filling factors. We find that the quasiholes are slightly 
better screened in the continuum than in the lattice and 
that they have the expected braiding properties. 

The paper is structured as follows. In Sec.jnj we intro¬ 
duce a family of wave functions and corresponding Hamil¬ 
tonians. The wave functions are bosonic and fermionic 
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Laughlin states on lattices for particular choices of the 
parameters. In Sec. |III[ we show that by choosing the 
parameters differently, it is possible within the same set 
of models to practically reproduce the continuous limit 
of the Laughlin states. This gives us a strategy to in¬ 
terpolate between the lattice and continuum Laughlin 
states in a way, where we know the Hamiltonian as well. 
In Sec. IV, we describe how we do the interpolation in 
practice, and we provide numerical results demonstrat¬ 
ing the applicability of our method. The models also 
allow Laughlin quasiholes to be added, and we investi¬ 
gate their properties in Secs. [V| and VI In particular, we 
show that the size of the quasiholes converges when going 
towards the continuum limit and that the quasiholes are 
slightly larger in the lattice than in the continuum. We 
also argue that the quasiholes have the expected braiding 
properties. In Sec. Ivnl we compute the effective mag¬ 
netic field seen by the quasiholes numerically. Finally, 
we conclude the paper with a discussion of the results in 
Sec. lyml 


II. WAVE FUNCTION AND HAMILTONIAN 

In this section, we introduce the lattice model that we 
are considering. We first define a two-dimensional lattice 
with sites, in the complex plane, at the positions = 
1,..., Vq. Laughlin-like states can be constructecP^EUon 
this lattice in terms of the vertex operators, Vn {zj) = 
(_1)0-I)"j ; exp[i{qnj - l)(l){zj)/^] with ; : 

the standard normal ordering in Conformal Field The¬ 
ory (CFT]P3, (p^z) a massless chiral free boson, q a 
positive integer and rij G {0,1} the number of hard¬ 
core bosons/fermions at lattice site Zj for q even/odd. 
Quasihole excitations over the Laughlin states can be 
made at arbitrary complex positions rji, i = 1,...,(5, 
by means of the quasihole vertex operators, Wp.{r]j) = 
: ex.-p[ipj(l){r]j)/-y/q] :, where pj is a positive integer that 
determines the charge of the quasihole as will be clear in 
a moment. The lattice wave functions are given by 

1 ^ 9 )= '^q{m^Q,ni^No)\ni,---,nNo), i'^) 

ni,...,nNQ 

where tji^q is shorthand notation for 771 ,..., ? 7 g, 

^qiVi^Q, ni^No)oc (0| Wp, (7?i) Wp 2 ( 772 )... WpQ (tiq) 

Vi (^ 1 ) Kx2 (^ 2 ) • • • i^No) |0) (2) 

and |0) is the CFT vacuum. Eq. ([^ evaluates tcP 

C{vi^Q)-^5n l[{zi - z,r^^^ llivi - 

i<j i<j 

X (3) 

with C a real normalization constant and = 1 for 
rii = {No - Yl?=iPi)h and (5„ = 0 otherwise. The 


lattice filling fraction is defined as z^Lat = ^i)/^o and 

in the absence of quasiholes {pi = 0 for all i) it becomes 
the Landau level filling fraction, z/Lat = ^ = 1/^- Note 
that from the definition of z/Lat and using the previous 
constraint on rii the charge of the j-th quasihole is 
given by Qj = Pj/q. 

If the lattice is defined on a disc-shaped region and 
the area per lattice site, a, is the same for all sites, 
then it can be showrP^ that the norm of the fac¬ 
tor's and Y{i^j{zi - Zj)~^^ approaches 

the usual Gaussian factors exp( —and 

exp(—in the thermodynamic limit. 
The phase factors can be transformed away if desired, 
and the state ^ then looks very similar to the Laughlin 
state in the continuum limit 


^qim^Q, ni^No) = C'im^Q) 

i<j 


where we have absorbed in C' the products — 

rjr^.yiPjlQ and fli ~ that are independent of 

Ui. Indeed, the state @ is a lattice version of the Laugh¬ 
lin state in the continuumi.e. in the absence of 
quasiholes it has a uniform lattice density and its quasi¬ 
hole excitations have the same topological properties as 
the Laughlin quasiholes in the continuum. Note that 
from the exponentials in Q we can define the lattice 
magnetic length, ^ a 

Using the results in [12] and |T6l it is possible to con¬ 
struct a Hamiltonian 


No 

Hq = Y,aIn 

i=l 


-I 2 


2 VU -j 

^ni- -(ATo-^p,) 

i = l ^ j = l 

Aj = ^ ——- [dj - di{qnj - 1)] - ^ p di 


with 




i=i 


Vj 


( 5 ) 


whose ground state is Here, c is a positive constant, 
and we note that the second term in Hq simply fixes 
the number of particles in the ground state to — 

{No — Ylf=i Pj)/Q' We find numerically for small systems 
that the ground state of Hq is unique for random choices 
of the lattice and quasihole coordinates. 


III. THE CONTINUUM LIMIT 

In order to approach the continuum limit, we need the 
number of lattice sites per area to increase, while keeping 
the number of particles per area constant. In [T^ and [T3l 
it has been demonstrated that this can be achieved by 
replacing the vertex operators Vnj{zj) by the vertex op¬ 
erators V^.{zj) = : exp['z(gnj — C)^(^i)/v^] ^ 
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where ^ is a positive number. Increasing the number of 
lattice sites to = A^o + AA^, this modifies the neutral- 

ity condition to “ J2?=iPvi)/‘l- We can 

hence keep the number of particles constant by choosing 
^ = Nq/N, and the continuum is obtained for N ^ oo. 

Here, we show that we can obtain a continuum limit 
of the lattice model (§ and by adding to it AN 
No uniformly distributed extra lattice sites and the same 
number of uniformly distributed background quasiholes 
with charge 1/q without changing the area of the system. 
In this way, the number of particles rii = {Nq — 

Ylf=iPr]i)/Q again stays the same. The advantage of this 
approach compared to the other is that the Hamiltonian 
^ remains valid. Note that the background quasiholes 
are placed at fixed positions with the only scope to reach 
the continuum limit, and they therefore play a different 
role than the physical quasiholes. To easily distinguish 
the two, we shall use pj to denote the positions of the Q 
physical quasiholes and Wj to denote the positions of the 
AN background quasiholes. 

We first discuss the continuum limit in the absence of 
physical quasiholes, i.e. Q = 0. From Eqs. ^ and ([^, 
in the presence of AN background quasiholes and AN 
extra sites, the lattice Hamiltonian and its ground state 
become 


N 


Han = 


i=l 


N 


- No/q 


_i=l 


with 


N ^ AN ^ 

Ai= N - [dj - di{qnj - 1 )] - Y] - d. 

^ ^ ^^ — ID ■ 




i=i 


(6) 


and 


Ylizi - - Wj)^ 

i<j i<j 

hi hi 

Note that the factors -{wi—Zj)^^ and 
can be written as 

hi 

Ei rij In \wi-Zj I y];. '>^i in \zi-Zj I ^icx ^ 

with a a phase factor that can be transformed away from 
the wave function if desired. In the following, we assume 
that both the lattice sites and the background quasiholes 
are distributed uniformly. We shall use the term unit 
cell to refer to the unit cell of the original lattice with 
AA/" = 0, and we shall take a to mean the area of this 
unit cell. In the limit AN oc, the sums in the last 
member of Q can be written as J]] • Uj \n\wi — Zj \ 

Pw Yjj rij f dw In \w - Zj \ and Y^j^i rii \n\zi - zj \ 


PzYil dzniln\zi - z\ with = n^/a and = n^/a 
being the densities of sites and quasiholes and Uz and 
the corresponding number of sites and quasiholes within 
one unit cell. Assuming in addition that the lattice is 
defined in a disc-shaped region, Eq. <§ turns into, 

Yliwi - ZjN Ylizi - = eP" Y-i 

hj iT^j 

^-Pz Eili s dzm In \zi-z\^ia _ ^2Tv{py,-p^) 

( 9 ) 

^Erom ^ it is clear that in order to obtain the right 
Gaussian factors (see Eq. @) in the continuum limit we 
must require that + 1/a, i.e. 

Tlz — T !• 

This equation states that the number of background 
quasiholes must be equal to the number of extra lattice 
sites. Einally, because the factors and 

Yli ~ independent of Ui they can be ab¬ 

sorbed into the normalization constant C and we find 
that the state 0 can be written as the state Q, with 
Pi = 0 for all i, but with the lattice sites covering the 
full space. Therefore, it is a continuum version of the 
lattice Laughlin state §-@. Erom here on we will use 
the terms lattice and continuum limit for AN = 0 and 
AN /Nq 1 respectively. 

The Laughlin states in the continuum are known to be 
the exact and unique ground state of the Hamiltonian 

H = Ho + V<-'^\ (11) 

with Ho the Landau level Hamiltonian and with 
the Trugman-Kivelson short-range potential interaction 
given by,^ 


N 

V(i) = ^ _ zj) (12) 

i<j 

for g' odd and 

= ( 13 ) 

i<j 

for q even. However it is important to remark that our 
Hamiltonian 0 is non-local and therefore very different 
in nature from the local one in Eq. 0. 

Einally, we can construct in a similar way a set of physi¬ 
cal quasiholes pi,..., pg in the continuum limit described 
by the ground state 

^q,AN,Q ~ ^'g.AiV P(M^i “ Vj)^ [[{Vi “ Vj)~ 
hi i<j 

X Yiivi - zjY''^i 

hi hi 
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of the Hamiltonian 


N 


HaN,Q — ^ 

i=l 

^ 1 

= ^2 ~ 


N 




A'o-ELft 


Li=l 


i=i 




with 


(15) 


and Ki defined in 


a) 



subcell 



IV. CONSTRUCTION TO OBTAIN A 
ROUGHLY UNIFORM PARTICLE DENSITY 


As shown in the previous section, to reach the con¬ 
tinuum limit we have to add to the original lattice an 
additional set of AN ^ Nq uniformly distributed lattice 
sites and the same number of uniformly distributed quasi¬ 
holes so that (10) is fulfilled. There are of course many 


different ways to reach such a limit. Naively one could 
consider, for instance, a random distribution of lattice 
sites and quasiholes satisfying (10). However, because 


from a computational point of view it is not possible to 
have a perfectly uniform distribution, many of the lat¬ 
tice sites will be close to the quasihole positions and the 
probability to find a particle on those lattice sites will 
be very small. This means that the lattice density per 
unit cell will be nonuniform and indeed we have checked 
that it is nonvanishing only in a small region of the cells. 

To overcome this problem and in order to interpolate 


.... H 
+ + + + + + '-' 


I- + 

+ 7 # + 

+ -I- 


FIG. 1: a) Lattice limit (AAc = 0) and unit cell, b) Poor 
continuum limit approximation (ANc = 7): unit cell made 
of 8 sites and 7 quasiholes (as indicated inside the unit cell). 


between the lattice and the continuum limit we propose 
here to consider the lattice structure in Fig. with 8 
sites distributed along the sides of the unit cell and 7 
quasiholes all of them placed at the center of the unit 
cell (note that Uz = ^ 1)- This is clearly a poor ap¬ 

proximation of the continuum limit because the number 
of extra lattice sites in a single unti cell, A Ac, is very 
small (AAc = 7). However, this construction can now 
be used as a building block to reach the continuum limit, 
as shown in Fig. where now the unit cell is made of 
3x3 subcells, each of them with the same structure of 
the unit cell of Fig. [^. Note that in order to satisfy the 
constraint (10) we keep the same number of quasiholes 


and lattice sites in all the subcells, except in the central 
one, that contains 7 quasiholes and 8 sites as shown in 


FIG. 2: Continuum limit approximation (AAc = 71) : a) 
Unit cell made of 9 subcells with the same structure as the 
unit cell in Fig. All the subcells contain 8 quasiholes 
except for the one in the center with 7 quasiholes (note that 
this setting satisfies the constraint Uz = riw 1)- The di, d2 
and da parameters are tuned by a Metropolis algorithm in 
order to reach an approximately uniform density in every 
subcell, b) Full lattice made of the unit cell in a). 


Fig. 1^. We also introduce, in every unit cell, a set of 
parameters di, d 2 and ds (see Fig. i^) to be tuned by a 
Metropolis algorithm, in order to better screen the effect 
of having one quasihole less in the subcell in the center 
than in the other subcells and thereby reach a more uni¬ 
form density. In Fig. we show, for the 3x3 sub cell 




Distance along a: axis 



Distance along x axis 


FIG. 3: a) Two-dimensional density as defined in eq.( |16| ) for 
the 3x3 subcell case of Fig. [^(AAc = 71) and for u = 1/3. 

b) Subcell density along the x-axis of the plot in a). The 
density is almost uniform with small oscillations around the 
mean value z//9 ~ 0.0370. The tuning parameters are given 
(in a=l unities) by: di = 2/9 + 0.041087, 
d 2 = 4/9 - 0.011062 and ds = 8/9 - 0.001649. c-d) Same 
plots as in a)-b) but for the lattice limit. 


case and u = 1/3, the two-dimensional subcell density 
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defined as, 


l{Xa,y/3)= ^ (rii), 

ieixc,,y^) 


(16) 


with {xa^y^) a set of coordinates pointing to the sub¬ 
cells and where the sum runs over all lattice sites i within 


the {xa^yp) subcell (note that in the lattice limit eq.(16) 
becomes the usual lattice density). Observe that the 
density is very close to uniform. This is more clear 
from Fig. ^jp where the subcell density along the x-axis 
presents small oscillations around its mean value u/9 (the 
sub cell density is given by u divided by the number of 
subcells per unit cell). In Fig. [^-d we plot the same 
densities but for the lattice limit obtaining the expected 
results.l^ As we will show later a good approximation of 
the continuum limit is given by a unit cell made of 5 x 5 
subcells, i.e. a unit cell with ANc = 199 sites and 5 
tuning parameters, chosen in a similar way as in Fig. [^. 
The density for this case and for u = 1/3 is shown in 
Fig. El with the mean density given by u/2b. Note that 
the density oscillations (Fig. El^) are higher than for the 
3x3 subcell case (Fig. |^) because we have to screen 
a much larger number of quasiholes with a few tuning 
parameters. This could be improved by including more 
parameters but, as it is clear from Fig. with this min¬ 
imal setting the oscillations are relatively small and the 
density of the system is close enough to the uniform case. 



Distance along 3) axis 



FIG. 4: a) Two-dimensional density as defined in eq.|T^ for 
the 5x5 subcell case {ANc = 199) and for u = 1/3. b) 
Subcell density along the x-axis of the plot in a). The 
density is close to uniform with small oscillations around the 
mean value iy/2b ~ 0.0133. The tuning parameters are given 
(in a=l unities) by: di = 2/15 — 0.012522, 
d2 = 4/15 - 0.011956, ds = 8/15 - 0.002250, 
d 4 = 10/15 - 0.000954 and ds = 14/15 + 0.000027. 


V. QUASIHOLE RADIUS AND CHARGE 

To study the extent of the quasiholes and their charge 
we split one particle into two quasiholes of charges 1/q 
and {q — l)/q and we place them far away from each 
other. To verify that the quasihole charge is equal to 1/q 


we compute the excess charg^^SHU 


^ 1 
[/’qh(ri) - p{ri)] ; ri = (i + -)e, i = 0,1, 2,... 

i=0 

(17) 

with pqhi'^i) the radial density at distance of one of 
the quasiholes, 

Pqh{ri)= (18) 


(n/)qh the quasihole mean density at position 2 )/, e a small 
real number, 77 the quasihole position and with p(r^) in 
Eq. ( prl ) given by (18) but replacing (^0* 



2D quasihole density 



Radial distance from quasihole 





FIG. 5: a,c,e): Quasihole subcell density, 
Tiqh{xa,xp) — n{xa,xp), for the lattice limit (a) and the 
3 X 3 (c) and 5 x 5 (e) subcell cases and for u = 1/3 
quasiholes with charges 1/3 and 2/3 placed at positions 
771 = —4 and 772 = 4, respectively. b,d,f): Excess charge, Q, 
of both quasiholes, for v —l/Z and again for the lattice limit 
(b) and the 3 x 3 (d) and 5 x 5 (f) subcell cases. Note that 
Q tends to the quasihole charges 1/3 (red curve) and 2/3 
(blue curve), as expected. 


In Eig. |^,c,e we show, for z/ = 1/3 and for the lattice 
limit and the 3x3 and 5x5 cases, the sub cell density, 
nqh(a;a,a;/ 3 ) - n{xa, xp) (see Eq. <E1I)> with nqh{xc„xp) 
the subcell density including two quasiholes with charges 
1/3 and 2/3 at positions 771 = —4 and 772 = 4, respec¬ 
tively. Erom the figures it is clear that there is a very 
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Lattice {ANc — 0) 

3x3 (AA^e = 71) 

5x5 (AAe = 199) 

3.0(3) 

2.0(3) 

2.0(5) 


TABLE I: Quasihole radius, i^qh, for quasiholes with 
charges l/3atz/=l/3 and for different values of ANc. 


good screening of quasiholes for different values of ANc. 
In Fig. ,d,f we place the quasiholes far away from each 
other and we show, again for z/ = 1/3 and for the previous 
cases, that the excess charge, Q, tends to the expected 
quasihole charges 1/3 and 2/3, respectively. 

To compute the q uasih ole radius, i?qh, we use the sec¬ 
ond moment of 


-^qh 


E£rlpqh(ri)-p(r,)|rf 

T;iZoMn)-pin)\ ’ 


(19) 


with jmax defined in such a way that is a distance far 
away from any quasihole, where Pqh(^jmax) “P(fomax) 
ishes. In table [Ipl| we show the values of i?qh in magnetic 
length units, r^for the quasiholes introduced previously 
with charges 1/3 and 2/3 at u = 1/3 and for the lattice 
limit and the 3x3 and 5x5 subcell cases. Observe that 
i?qh decreases with ANc because of the better screening 
when one approaches the continuum limit (see Fig. 
and that the values of i?qh for the 3x3 and 5x5 cases 
are very close, indicating that the 5x5 subcell case is a 
good approximation for the continuum limit. 

Finally we want to remark that the values we obtain for 
i?qh in the continuum limit are v ery clo se to the ones ob¬ 
tained in the standard FQH effeclP^^ and other systems 
containing Laughlin-like quasihole states like Fractional 
Chern Insulators.!^ 


Lattice {ANc — 0) 

3x3 (AAc = 71) 

5x5 (AAc = 199) 

3.0(1) 

2.3(5) 

2.3(5) 


TABLE II: Quasihole radius, i^qh, for quasiholes with 
charges 2/3atz/=l/3 and for different values of ANc. 


VI. QUASIHOLE BRAIDING 


In this section we show that the anyonic statisticJ^^^^ 
remains invariant in the interpolation between the lattice 
and the continuum limit. 

To compute the anyonic statistics we proceed as in the 
previous section and we split up one of the particles into 
two quasiholes 771 and 772 , with charges 1/q (p^^ = 1) 
and {q — l)/q {Pr ]2 = ^ ~ !)• The anyon ic statistics, 7 , is 
defined in terms of the Berry phasd^^^ and monodrom}!^ 
(i.e. the change obtained from analytical continuation of 
(14) when the quasiholes move around) as. 


27r7 = + arg{Ma) - 0^ - arg{Mb), (20) 


with Oa and Mq the Berry phase and monodromy due 
to the braiding of the quasiholes (Fig. |^) and 0^ and 
Ml) the same quantities but without quasihole braid¬ 
ing (Fig. 1^). The Berry phase for the state (14) can 



FIG. 6: Set up for anyonic statistic computation of 
quasiholes 771 and 772 with charges t and In a) the 
quasihole 772 is placed inside the curve L (at position 772 ) cind 
in b) far away from it (at position 772 ). 


be computed directly from the normalization constant 
C{7]i,7]2,Wi^An)^ 



1 dC^ 

drii 


drii - 



1 dC^ 


drii 


( 21 ) 


and using the explicit form of the wavefunction it can be 
written as 


0 = —Im 


1 f drii ^ 

q ^ Jr Vi -Wi q Jr Vi- m 

1 7 dljl ^ f {rij) 

q , Jr Vi - Zi , Jr 01- Zi 


dr]i 


( 22 ) 


Observe that the first three members involve analyti¬ 
cal functions and they can be computed easily. Using 
Eq. (p2|, 


0a-db = i\ 27ri 


.q-1 




{ni)a - {ni)b 
01 - Zi 


drti 


(23) 


where {ui)^ ^he density on the i-th lattice site 

with the quasihole at position 772 ( 72)5 shown in Fig.lm 
Note that from the numerical results of Sec.V (see Fig.|^ 
the function {rii)^ — {rii)^ in (23) is zero at every lattice 
site except around the positions 772 and 772 , i.e. 


(«i)a - {^i)b = -r{Zi - 02 ) + f{Zi - T]^), (24) 


with /" (/^) some analytical function, nonvanishing only 

around 772 ( 72)5 ^^^h that ~ = {q~ 

l)/q. Finally, from Eq. ( [^ and Eq. ( [^ it is clear that 
Oa — Ob = 0 and therefore we obtain from Eq. ( [2Q| ) that 
7 = {arg{Ma) — arg{Mb))/2 t: = {q — l)/q and thus the 
quasihole anyonic statistics can be read directly from the 
monodromy and it remains invariant along the lattice- 
continuum limit interpolation. 
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VII. MAGNETIC FIELD 


The Berry phase acquired when a quasihole moves 
around a closed loop while all other quasiholes are far 
away can be interpreted as an Aharonov-Bohm phase of 
a charged particle in an effective magnetic field i.e. 


0 = -iln(M) 


{pie/q){2TTh)a 



B{xi,yi)dxidyi. 


(25) 


In the lattice limit the effective magnetic field felt W 
a quasihole in a unit cell is not completely uniforrrP^. 
However, in this section we show that it approaches the 
uniform value B ^ —h/ (ea) (note that the magnetic flux 
through a unit cell is —hje) as we approach the contin¬ 
uum limit. Using ( [2^ and the divergence theorem on 
@ the magnetic field, H, can be related to the occupa¬ 
tion number (in natural units h = e = 1 and taking 
a = 1), 


B = - 


qPm {riiTij) - {fij) {rij) ^ 
TT ^ {r]a - Zi){f]a - Zj)' 


a = 1,2 (26) 


where again we split up a particle in two quasiholes at po¬ 
sition 771 , 772 , far away from each other, and with charges 
1 /^ (Pm = 1 ) and (^ - 1 )/^ {pr^^ = ^ - 1 ) respectively. 
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FIG. 7: a) Effective magnetic field, B, felt by a z/ = 1/3 
quasihole with charge 1/3 in a 10 x 10 grid (as explained in 
the main text) of a single unit cell, for the 3x3 subcell case, 
b) One-dimensional version of a) for the dependence of B in 
every subcell grid. Note that the magnetic field is almost 
uniform with small oscillations around the uniform value -1 
(in units h = e = c = 1 and a = 1). 


In Fig. we show the dependence of the magnetic 
field felt by a = 1/3 quasihole with charge 1/3 for the 
3x3 case. We divide a unit cell in a 10 x 10 grid, and 
in Fig. we plot the two-dimensional dependence of 


B on the grid. Fig. is a one-dimensional version of 
Fig.[7fi and shows the dependence of B in every subcell 
of the grid to make more clear that the magnetic field 
is almost uniform with small physical fluctuations (we 
have checked that the Monte Carlo errors are smaller 
than the fluctuations) of order 0.001 around the mean 
value —1. Finally, we have checked that similar results 
hold for the quasihole with charge 2/3 and for both 
quasiholes in the 5x5 subcell case, though in this case 
the Monte Carlo errors are bigger and it is more difficult 
to distinguish the physical magnetic field fluctuations 
from the error bars. 


VIII. CONCLUSION 

In conclusion, we have proposed how one can use CFT 
to construct both lattice and continuum models with 
Laughlin-like ground states and to interpolate between 
the two limits. We have also shown that the topological 
properties of the models are as expected and computed 
the size of the quasiholes. 

An interesting feature of the models is that both the 
Hamiltonian and the unique ground state are known an¬ 
alytically, both with and without quasiholes. This allows 
us to compute a number of relevant properties easily with 
Monte Carlo simulations. It also allows us to show an¬ 
alytically that the braiding properties of the quasiholes 
are as expected if the quasiholes are screened. The mod¬ 
els hence remain within the same topological phase for 
all situations, where screening occurs. 

The models are also interesting, because they allow us 
to compare the Hamiltonian obtained from CFT in the 
continuum limit with the usual delta function interaction 
Hamiltonian for the Laughlin states. The Hamiltonians 
turn out to be different, so that the CFT approach pro¬ 
vides a different set of models displaying FQH properties. 

Finally, we note that the approach presented in this 
paper can also be used to interpolate other FQH states 
constructed from CFT between the lattice and the con¬ 
tinuum limit. This is so because the background charge 
always appears in the same way in the CFT correlators 
used to construct the states. 
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